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Abstract 
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1. Introduction 



Forty years after Pontecorvo's discovery of the theory of neutrino oscilla- 
tions IH, there are today three experimental indications in favor of the existence 
of this phenomenon. The strongest evidence is believed to come from the at- 
mospheric neutrino anomaly observed by the Super-Kamiokande experiment [Q, 
following the previous measurements of the Kamiokande, 1MB and Soudan at- 
mospheric neutrino experiments ||. The atmospheric neutrino anomaly consists 
in a suppression of the ratio of measured atmospheric and v e with respect to 
the calculated one (see, for example, Ref. [||]). Furthermore, the new data of the 
Super-Kamiokande experiment || indicate a variation of the atmospheric z/ M flux 
as a function of the zenith angle, which correspond to a variation as a function 
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of the distance between the region in the atmosphere where the neutrino is pro- 
duced and the Super-Kamiokande detector where the neutrino is detected. This 
is a strong evidence in favor of neutrino oscillations and seems to be supported 
by the measurement of the flux of neutrino-induced upward-going muons of the 
MACRO experiment 0. 

The oldest indication in favor of neutrino oscillations comes from the so- 
lar neutrino problem (see, for example, Ref. ||), i.e. the suppression of the flux 
of solar i/ e 's observed in five solar neutrino experiments (Homestake, Kamiokande, 
GALLEX, SAGE and Super-Kamiokande 0) with respect to the calculated one ||. 

Finally, the LSND accelerator experiment has reported evidence in favor 
of — > v e and — > v e oscillations |J. 

Neutrino oscillations can occur only if neutrinos are massive particles, if 
their masses are different and if neutrino mixing is realized in nature. In this 
case, the left-handed flavor neutrino fields v aL [a = e, /i, r) are superpositions 
of the left-handed components u^l (k = 1, • • • , n) of the fields of neutrinos with 
definite masses m^. v aL = Y^k=i ^ ak VkL > where U is a unitary mixing matrix. 
The general expression for the probability of v a — > vp transitions in vacuum is 



p 



k=l ^ 



(1) 



where Am 2 ^ = — mj, L is the distance between the neutrino source and 
detector and p is the neutrino momentum. 

2. Four neutrino mixing and oscillations 

The three experimental indications in favor of neutrino oscillations re- 
quire the existence of three different scales of neutrino mass-squared differences: 
Am s 2 un ~ l(r 5 eV 2 (MSW) or Am 2 sun ~ 10~ 10 eV 2 (vacuum oscillations), Am 2 tm ~ 
5 x 10 _3 eV 2 and Am| BL ~ leV 2 , where Am| BL is the neutrino mass-squared 
difference relevant for short-baseline (SBL) experiments, whose allowed range is 
determined by the positive result of the LSND experiment. 

Three independent mass-squared differences require at least four massive 
neutrinos. The number of active light flavor neutrinos is known to be three 
from the measurement of the invisible width of the Z-boson, but there is no 
experimental upper bound for the number of massive neutrinos (the lower bound 
is three). In the following we consider the simplest possibility of existence of four 
massive neutrinos [ ID| , |TT, [1^, |1^, |14], [Uj. In this case, in the flavor basis, 
besides the three light flavor neutrinos v e , v T that contribute to the invisible 
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width of the Z-boson, there is a light sterile neutrino u s that is a SU(2) L singlet 
and does not take part in standard weak interactions. 

Two years ago we have shown [HJ that among all the possible four-neutrino 
mass spectra only two are compatible with the results of all neutrino oscillation 
experiments: 
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and (B) m x < m 2 <C m 3 < m 4 . (2) 
s / 
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In these two schemes the four neutrino masses are divided in two pairs of close 
masses separated by a gap of about 1 eV, which provides the mass-squared differ- 

m\ — m\ that is relevant for the oscillations observed in 
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the LSND experiment. In scheme A Am 
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explanation of the atmospheric neutrino anomaly and Amg Un = Am^ = m\ — m\ 
is relevant for the suppression of solar i/ e 's, whereas in scheme B Am 2 tm = Am| 3 



and Am? 



Am^. 



Let us define the quantities c a , with a = e, fi, r, s, in the schemes A and B 



as 
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\U a k\ 



(B) 



Cry 



\u a k\ 



(3) 



fc=l,2 



fc=3,4 



Physically c a quantify the mixing of the flavor neutrino v a with the two massive 
neutrinos whose Am 2 is relevant for the oscillations of atmospheric neutrinos (z/ 1; 
V2 in scheme A and z/ 3 , z/4 in scheme B). 

The probability of v a — > vp transitions (/3 7^ a) and the survival probability 
of v a in SBL experiments are given by [12] 
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with the oscillation amplitudes 
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where the index k runs over the values 1,2 or 3,4. The probabilities (^) have the 
same form as the corresponding probabilities in the case of two-neutrino mixing, 
P Va ^ Vp = sin 2 (20) sin 2 (Am 2 L/4p) and P Va ^ Ua = 1 - sin 2 (20) sin 2 (Am 2 L/4p), 
which have been used by all experimental collaborations for the analysis of the 
data in order to get information on the parameters sin 2 (29) and Am 2 (9 and Am 2 
are, respectively, the mixing angle and the mass-squared difference in the case of 
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two- neutrino mixing). Therefore, we can use the results of their analyses in order 
to get information on the corresponding parameters A a .p, B a . a and Am| BL . 

The results of all neutrino oscillation experiments are compatible with the 
schemes A and B only if [[13] 

c e < a° e and c M > 1 — a° , (6) 

where 



1 - J I - B° (a = e, fi) (7) 



ol 2 V V j 

and B®. a is the upper bound for the amplitude of v a — >v a oscillations obtained 
from the exclusion plots of SBL reactor and accelerator disappearance experi- 
ments. Hence, the quantities a° e and a° depend on Am 2 . The exclusion curves 



obtained in the Bugey reactor experiment [17] and in the CDHS and CCFR ac- 



celerator experiments [0 imply that both a° e and a are small [pTQf] : aP e < 4 x 10 2 



and < 2 x 10 1 for any value of Am 2 in the range 0.3 < Am 2 < 10 3 eV 2 . 



The smallness of c e in both schemes A and B is a consequence of the 



solar neutrino problem [p^| . It implies that the electron neutrino has a small 
mixing with the neutrinos whose mass-squared difference is responsible for the 
oscillations of atmospheric neutrinos (z/i, i>2 in scheme A and v^, 1/4 in scheme B). 
Therefore, the transition probability of electron neutrinos and antineutrinos into 
other states in atmospheric and long-baseline (LBL) experiments is suppressed. 
Indeed, it can be shown [TJJ that the transition probabilities of electron neutrinos 
and antineutrinos into all other states and the probability of u^— >v e transitions 
in vacuum are bounded by 

1 - P^] < al (2 - a® (8) 
^ B ( L J ) <minfa°(2-a°) ,a° + i< ;e ) (9) 

where A°^. e is the upper bound for the amplitude of v^— >v e transitions measured 
in SBL experiments with accelerator neutrinos. 

The upper bound (||) obtained from the 90% CL exclusion plot of the 
Bugey experiment |17| is shown in Fig. [I] (solid line). The shadowed region in 
Fig. HI corresponds to the range of Am 2 allowed at 90% CL by the results of the 
LSND and all the other SBL experiments. The upper bound for the transition 
probability of ^ e 's into other states obtained from the recently published exclusion 
plot of the CHOOZ |2D] experiment is shown in Fig. [I] (dash-dotted line). One 



can see that the CHOOZ upper bound agrees with that obtained from Eq.(^). In 
Fig. [I] we have also drawn the curve corresponding to the expected final sensitivity 
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Figure [I] Figure ^| 

of the CHOOZ experiment (dash-dot-dotted line). One can see that the LSND- 
allowed region imply an upper bound for 1 — P§^^ e of about 5 x 10~ 2 , which is 
smaller than the expected final sensitivity of the CHOOZ experiment. 

The bound obtained from the 90% CL exclusion plots of the Bugey 



is 



experiment and of the BNL E734, BNL E776 and CCFR experiments [pi 
depicted by the dashed line in Fig. §. The solid line in Fig. || shows the upper 
bound on F^f'l taking into account matter effects. The expected sensitivities of 
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the K2K long-baseline accelerator neutrino experiment is indicated in Fig. g 
by the dash-dotted vertical line. The shadowed region in Fig. || corresponds to 
the range of Am 2 allowed at 90% CL by the results of the LSND and all the other 
SBL experiments. It can be seen that the results of SBL experiments indicate 
an upper bound for E_\ B }_s smaller than 4 x 10 -2 and smaller than the expected 

sensitivity of the K2K experiment. 

The two schemes A and B have identical implications for neutrino oscilla- 
tion experiments, but very different implications for neutrinoless double-/? decay 
experiments and for tritium /3-decay experiments. Indeed, in scheme A 
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whereas in scheme B 
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where (m) = rrii is the effective Majorana mass that determines the 

matrix element of neutrinoless double-/? decay and m( 3 H) is the neutrino mass 
measured in tritium /3-decay experiments. Therefore, in scheme B |(m)| and 
m( 3 H) are smaller than the expected sensitivity of the next generation of neu- 
trinoless double-/? decay and tritium /3-decay experiments. The observation of a 
positive signal in these experiments would be an indication in favor of scheme A. 

Summarizing, the results of neutrino oscillation experiments indicate that 
only the two four-neutrino schemes @) are allowed and the electron neutrino has 
a very small mixing with the two massive neutrinos that are responsible for the 
oscillations of atmospheric neutrinos (z/i, V2 in scheme A and z/3, z/4 in scheme B). 
Hence, the two schemes have the form shown in Fig. 0, where v e is associated 
with the two massive neutrinos neutrinos that are responsible for the oscillations 
of solar neutrinos (z/3, z/4 in scheme A and z/i, z/2 in scheme B), with which it has 
a large mixing, whereas is associated with the two massive neutrinos neutrinos 
that are responsible for the oscillations of atmospheric neutrinos, with which the 
muon neutrino has a large mixing. The results of neutrino oscillation experiments 
do not provide yet an indication of where v T and v a have to be placed in the two 
schemes represented in Fig. |3|. In the following Section we will show that the 
standard Big-Bang Nucleosynthesis constraint on the number of light neutrinos 
provide a stringent limit on the mixing of the sterile neutrino with the two massive 
neutrinos that are responsible for the oscillations of atmospheric neutrinos [O, W 



3. BBN constraints on 4-neutrino mixing 

It is well known that the observed abundance of primordial light elements 
is predicted with an impressive degree of accuracy by the standard model of Big- 
Bang Nucleosynthesis if the number N v of light neutrinos (with mass much smaller 
than 1 MeV) in equilibrium at the neutrino decoupling temperature (T^ ec ~ 2 MeV 
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for v e and T^ ec — 4 MeV for u^, u T ) is not far from three ||23|, p4f . 

The value of N v is especially crucial for the primordial abundance of 4 He. 
This is due to the fact that N v determines the freeze-out temperature of the weak 
interaction processes e + + n ^=> p + u e , v e + n p + e~ and n p + e~ + z/ e 
that maintain protons and neutrons in equilibrium, i. e. the temperature at which 
the rate IV (T) ~ GpT 5 (Gp is the Fermi constant) of these weak interaction 
processes becomes smaller than the expansion rate of the universe 



rr (T ,R(T) /87T3 T 2 

(Mp is the Planck mass), where R(T) is the cosmic scale factor and g* = 5.5 + 
1.75N U for m e < T < m^. If iV^ = 3 the primordial mass fraction of 4 He, 
Yp = mass density of 4 He / total mass density, is Yp ~ 0.24, which agrees very 
well with the observed value |4| Fp bs = 0.238±0.002. Since Yp is very sensitive to 



variations of N u , it is clear that the observed value of Y P puts stringent constraints 
on the possible deviation of A^, from the Standard Model value N u = 3. 

In the four-neutrino schemes @ standard BBN gives a constraint on neu- 
trino mixing if the upper bound for N u is less than four. In this case the mixing 
of the sterile neutrino is severely constrained because otherwise neutrino oscilla- 
tions would bring the sterile neutrino in equilibrium before neutrino decoupling 



25|, |27|, |8|, |29], 0, |3T], [TT|], leading to N v = 4. In particular, we will show that 
standard BBN with N v < 4 implies that c s is extremely small. 

The amount of sterile neutrinos present at nucleosynthesis can be calcu- 
lated using the differential equation 

dn Ua 1 



£ (i^Jcoiir^ (l-nj, (13) 



dT 2HT 

where n Us is the number density of the sterile neutrino relative to the number 
density of an active neutrino in equilibrium and T Ua are the collision rates of the 
active neutrinos, including elastic and inelastic scattering [p9[| , T Ue = 4.0 G%T 5 
and IV = T Ut = 0.7 T Ue . The quantities (P Ua ->u s ) coil are the probabilities of 
v& — > v s transitions averaged over the collision time t co n = l/T Uoi . Hence, also n Ug 
has to be considered as a quantity averaged over the collision time. 



Equation (13) describes non-resonant and adiabatic resonant neutrino tran- 
sitions if t osc <^ ^coii ^ ^ex P - The condition t osc t co n means that neutrino 
oscillations have to be fast relative to the collision time. The characteristic ex- 
pansion time of the universe t exp is given by t exp = 1/H where H is the Hubble 
parameter H = R/R, which is related to the temperature T by H = —T/T ~ 
0.7 (T/1 MeV) 2 s" 1 (this value can be obtained from Eq.(|T2"D with m e < T < m M 
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3). The relation T Ue /H ~ 1.2 (T/l MeV) 3 shows that for temperatures 
larger than 2 MeV the collision time is always much smaller than the expansion 



and N v ~ 



time 25 



Since by definition N v is the effective number of massless neutrino species 
at Tdoc, in order to get a constraint on the mixing of sterile neutrinos we need 
to calculate the value of n Vs at Td ec produced by neutrino oscillations. With the 
initial condition n Ua {Ti) = (Tj ~ 100 MeV), the integration of Eq.(|T3"D gives f30 



dec I 



with 



Ti 



2HT 



(14) 



a=e,fi,r 



Imposing the upper bound n Vs (T& ec ) < 5N 
F < |ln(l -5N)\. 



N u — 3 one obtains the condition 



For the calculation of F the averaged transition probabilities (P Va -n, s ) 



coll 



must be evaluated and the effective potentials of neutrinos and antineutrinos due 



to coherent forward scattering in the primordial plasma |32 

rp4 



V e = -6.02 G F p 



M 2 W 



= V, V, 



fl,T 



iV and V a = . 



(15) 



must be taken into account (in the absence of a lepton asymmetry the effective po- 
tentials of neutrinos and antineutrinos are equal). Here p is the neutrino momen- 
tum, which we approximate with its temperature average (p) ~ 3.15T, Gf is the 
Fermi constant, Mw is the mass of the W boson and £ = cos 2 6^/(2 + cos 2 9w) — 
0.28, where 9w is the weak mixing angle. The propagation of neutrinos and 
antineutrinos is governed by the effective hamiltonian in the weak basis 
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U ] + diag 



V,ZV,£V,0 



(16) 
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It is convenient to subtract from Hy/ the constant term p + m 1 /2p + £V, which 
does not affect the relative evolution of the neutrino flavor states, in order to get 
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C/t+diag (l-OV,0,0,-eF ■ (17) 



From this expression it is clear that the relative evolution of the flavor neutrino 
states depends on the three mass-squared differences and not on the absolute 
scale of the neutrino masses. The effective hamiltonian in the mass basis is given 
byH' M = WH$yU: 



0, Am 2 ,!, Am^, Am^ 



+ E/+diag (l-£)V,0,0,-£V 



U. 



In the mass basis the mixing has been transferred from the mass term to the 
potential term. In order to calculate the evolution of the neutrino flavors it is 
necessary to parameterize the 4x4 neutrino mixing matrix U. However, since the 
second and third rows and the second and third columns of the diagonal potential 
matrix in Eq . fll8|) are equal to zero, it is clear that the values of the second and 
third rows of U, corresponding to and u T , are irrelevant and do not need to be 
parameterized. Furthermore, since c e is small in both schemes A and B, it does 
not have any effect on neutrino oscillations before BBN and the approximation 
c e = is allowed. Hence, the 4x4 neutrino mixing matrix can be partially 
parameterized as 
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in scheme A and 
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l sin if sin x — sin if cos x — cos if sin 9 cos if cos 9 i 

\ 



u 



sin 9 



(20) 



i — cos if sin 9 cos if cos 9 sin if sin x — sin if cos x } 



in scheme B, with < ip < ir/2. Notice that c s = sin 2 if in both schemes A 
and B. The dots in Eqs.(|19D and (|20| ) indicate the elements of the mixing matrix 
belonging to the and v T rows (U^ and U ri with i = 1, ... ,4), which do not 
need to be parameterized. In Eqs.(|l^) and ( p0|) we have parameterized only the 
elements of the mixing matrix belonging to the v e and v s lines (U e i and U S i with 
i = 1, . . . , 4) in terms of the three mixing angles 9, x, f- It is clear that this partial 
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parameterization of the mixing matrix (with the approximation U e i = U e 2 = 0) is 
much easier to manipulate than a complete parameterization, which would require 
the introduction of 6 mixing angles and 3 complex phases. 

Notice that no complex phase is needed for the partial parameterization 
of the mixing matrix in Eqs.(|19"D and (p0[), because the elements U e i and U S i 
with i = 1, ... ,4 can be chosen real. Indeed, if we consider, for example, the 
scheme A, the line U S i with i = 1,...,4 and the element ?7 e 4 can be chosen 
real because all observable transition probabilities are invariant under the phase 
transformation U a j — > e tXa U a j e %Vi , where x a and yj are arbitrary parameters. In 
order to show that also Z/ e3 can be chosen real, we multiply the unitarity relation 



with U* 3 U s3 . 



The imaginary part of the resulting relation gives 



Im[[/ e3 C/; 3 C/ e * 4 C/ S 4] 



lm[U el U* sl U* e3 U s3 ] + lm[U e2 U* s2 U* e3 U t 



,s3j 



(21) 



In the approximation U e \ = U e 2 = the right-hand part of Eq.([21|) vanishes. 
Therefore, since U S 3, U e 4, C/ s 4 have been chosen to be real, also U e z must be real. 
Since the mass-squared differences have a hierarchical structure, Am] 3 -C 



Amj! <^ Am 2 
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in scheme A and Am^ 



<^ Aml 3 Arnjj in scheme B, 



the 



effective hamiltonian H' M can be diagonalized approximately taking into account 
only one of the three Am 2 's for different ranges of the temperature T. Then, it 



can be shown that [[11], [16| the condition F < \ ln(l — 5N)\ gives the bound 
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< |ln(l-57V)| , 
(22) 



with gL = c s in scheme A and gL = 1 — c s in scheme B. 



Both terms in the left-hand side of Eq. (P2|) are positive and must be small 
if SN < 1. The SBL term, depending on Amg BL , is small if either c s is small 
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or large, but the atmospheric term, which depends on Am 2 tm , is small only if 
c s is small. Indeed, if c s is close to one we have (U^,!!^) ~ (cos x, sin x) m 
scheme A and (^3,^4) ~ (cos x, sin x) in scheme B. This means that, in order 
to accommodate the atmospheric neutrino anomaly, sin 2 2% cannot be small. This 
is in contradiction with the inequality ( p2[) and we conclude that the bound (22) 
implies that c s is small. 

Since c s is small only non-resonant transitions of active into sterile neutri- 
nos due to Amg BL are possible in scheme A, as illustrated in Fig. |] where we have 
plotted the effective squared masses (obtained from a numerical diagonalization 
of the hamiltonian (|T6D) as functions of T 6 [y e does not have resonant transitions 
into or v T because we have chosen c e = 0). Hence, the conditions for the 
validity of Eq. ([13]) are satisfied and the SBL term in Eq.(|22|) gives the bound 

-1/2 



C, < 1.1 X 10" 



leV 2 



ln(l-cWV)|. 



(23) 



On the other hand, since c s is small, a resonance occurs in scheme B at the tem- 
perature T res = 16(Am| BL /l eV 2 ) 1 / 6 ^ - 2c s \ 1 ' 6 MeV, as illustrated in Fig. §. The 
condition SN < 1 implies that this resonance must not be passed adiabatically. 
In this case the conditions for the validity of Eq.(^) are not fulfilled and the 
SBL term of Eq.(p2|) does not apply. Using an appropriate formula |27| for the 
calculation of the amount of sterile neutrinos produced at the resonance through 



non-adiabatic transitions one can show Jl6j that the BBN bound on c s in scheme 
B is given by 

c, < 1.1 x 10" 



( Am| BL \ 


-1/2 


V leV 2 ) 





ln(l-5iV)|. (24) 

Figure |6] shows the values of the bounds (^3|) and (^) obtained from the 
LSND § lower bound Am| BL > 0.27 eV 2 for 0.2 < 5N < 1. One can see that 
standard BBN implies that c s is extremely small. 
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4. Conclusions 

The BBN constraint on c s derived in the previous Section implies that u s 
is mainly mixed with the two massive neutrinos that contribute to solar neutrino 
oscillations (1/3 and 1/4 in scheme A and V\ and z/2 in scheme B) and the unitarity 
of the mixing matrix implies that v T is mainly mixed with the two massive neu- 
trinos that contribute to the oscillations of atmospheric neutrinos. Adding this 
information to the two schemes of neutrino mixing depicted in Fig. ||] we obtain 
the schemes shown in Fig. 

The two neutrino mixing schemes shown in Fig. || have the following 
testable implications for solar, atmospheric, long-baseline and short-baseline neu- 
trino oscillation experiments: 

• The solar neutrino problem is due to v e — > v s oscillations. This prediction 
will be checked by future solar neutrino experiments that can measure the 
ratio of neutral-current and charged- current events [ |33|| . 

• The atmospheric neutrino anomaly is due to z/ M — > v T oscillations. This 
prediction will be checked by LBL experiments. Furthermore, the absence 
of Up — > v s atmospheric neutrino oscillations may be checked in the future in 
the Super-Kamiokande atmospheric neutrino experiment p3 . 



— > v T and v e — > v s oscillations are strongly suppressed in SBL experiments. 
With the approximation c s ~ 0, for the amplitude of — > v T oscillations 
we have the upper bound A^ T < (c^) 2 , that is shown in Fig. ^ (solid curve) 
together with a recent exclusion curve obtained in the CHORUS experiment 
(dash-dotted curve) and the final sensitivity of the CHORUS and NOMAD 
experiments (dash-dot-dotted curve) |35| . 



If these prediction will be falsified by future experiments it could mean that some 
of the indications given by the results of neutrino oscillations experiments are 
wrong and neither of the two four neutrino schemes A and B is realized in nature, 
or that Big-Bang Nucleosynthesis occurs with a non-standard mechanism |36| . 



In conclusion, we would like to emphasize that if the analysis presented 
here is correct and one of the two four neutrino schemes depicted in Fig. § is 
realized in nature, at the zeroth-order in the expansion over the small quantities 
c e and c s the 4x4 neutrino mixing matrix has an extremely simple structure in 
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which the u e , v s and u^, v T sectors are decoupled: in scheme A 



U 



( 

cos 7 sin 7 
— sin 7 cos 7 



V 



o 



o 



cos# sin# \ 





- sin 9 cos 9 



(25) 



and in scheme B 



U 



( cos 9 sin 9 



\ — sin 9 cos 9 



\ 

cos 7 sin 7 
— sin 7 cos 7 











(26) 



where 9 and 7 are, respectively, the two-generation mixing angles relevant in solar 
and atmospheric neutrino oscillations. Therefore, the oscillations of solar and 
atmospheric neutrinos are independent and the two-generation analyses of solar 
and atmospheric neutrino oscillations yield correct information on the mixing of 
four-neutrinos. 

S.M.B. acknowledge the support of the "Sonderforschungsbereich 375-95 
fuer Astro- Teilchenphysik der Deutschen Forschungsgemeinschaf " . 
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